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Abstract 

We perform an updated phenomenological analysis of the Mikheyev-Smirnov- 
Wolfenstein (MSW) solutions of the solar neutrino problem, assuming os- 
cillations between two and three neutrino families. The analysis includes 
the total rates of the Homestake, SAGE, GALLEX, Kamiokande and Super- 
Kamiokande experiments, as well as the day-night asymmetry and the 18-bin 
energy spectrum of Super-Kamiokande. Solutions are found at several values 
of the #i3 mixing angle. Among the most interesting features, we find that 
solar neutrino data alone put the constraint #13 < 55°-59° at 95% C.L., and 
that a fraction of the MSW solutions extends at and beyond maximal (y\, 1^2) 
mixing (#12 > 7r/4), especially if the neutrino square mass splitting is in its 
lower range (m| — ml ~ 10~ 7 eV 2 ) and if #13 is nonzero. In particular, bi- 
maximal (or nearly bimaximal) mixing is possible for atmospheric and MSW 
solar neutrino oscillations within the stringent reactor bounds on #13. 

PACS number(s): 26.65.+t, 13.15.+g, 14.60.Pq, 91.35.-x 
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I. INTRODUCTION 



It is widely recognized that the combined sources of evidence for neutrino flavor tran- 
sitions coming from the solar neutrino problem [0J and from the atmospheric neutrino 
anomaly @ demand an approach in terms of three-flavor oscillations among massive neutri- 
nos (ui, V21 ^3) 0. The three-flavor v parameter space is then spanned by six variables: 

5m 2 = m\ — m\ , (la) 
m 2 = m\ — m\ , (lb) 

^ = #12 e [0,tt/2] , (lc) 

= ^13 G [0,vr/2] , (Id) 
V = 023e[O,7r/2], (le) 
5 = CP violation phase , (If) 

where the 6ij rotations are conventionally ordered as for the quark mixing matrix 

In the phenomenologically interesting limit \5m 2 \ <C \m 2 \, the two eigenstates closest in 
mass (^1,^2) drive solar neutrino oscillations, while the "lone" eigenstate z/ 3 drives atmo- 
spheric neutrino oscillations. In such a limit (see p,|J and refs. therein): 1) the phase 5 
becomes unobservable; ii) the atmospheric parameter space is spanned by (m 2 ,ip,cf)); and 
Hi) the solar neutrino parameter space is spanned by (5m 2 , u>, </>).[] 

Building on a previous work ||, we perform a thorough analysis of the available solar neu- 
trino data in the (5m 2 ,uj,<j)) variables, in the context of the Mikheyev-Smirnov-Wolfenstein 
(MSW) oscillation mechanism 0. There are several motivations to revisit, update and im- 
prove the analysis in [Q, the most important being the need to include the high-statistics 
Super-Kamiokande (SK) observations of neutrino events (total rate, ||, energy spectrum 
|2|f7|], and day- night difference 0,0). Indeed, the few post-SK papers on three-flavor MSW 
oscillations we are aware of [PHTT|| include neither the spectrum information nor day-night 
variations. This contrasts with the more familiar two-family MSW analyses, which have 
been regularly updated with state-of-the-art global fits at different detector lifetime days; 
see, e.g., the SK official fits ( [P| |i2"| , |i"5[ and refs. therein), as well as various analyses by 



independent research groups: |TJ| (100 days), |L5| (300 days), |L6| (500 days), and [|T7|JT1 
(800 days). 

After the work ||, there have been also other relevant experimental and theoretical 
improvements (implemented in the present analysis): updated measurements of the exper- 
imental (total) rates in the chlorine (CI) |H| and Kamiokande (K) |20| experiments and in 



the gallium (Ga) detectors SAGE [glj and GALLEX |22]; new standard solar model (SSM) 
estimates for the neutrino fluxes and their uncertainties, and for the neutrino production 
regions in the Sun [p3|j ; updated calculations of the 8 B neutrino spectrum |24|], of the u, 



scattering cross section on electrons [J25|, and of the u e neutrino absorption cross section 



1 In the special case (f) = 0, the atmospheric and solar parameter spaces are decoupled into the 
two-family oscillation spaces (5m 2 , uj) and (m 2 , ■*/>). 
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in chlorine |24]] and in gallium [26[|; analytical methods to compute the (time-averaged) v e 
survival probability in the Earth . 



In addition, some relatively old topics, such as the role played by a possibly large hep 
neutrino flux p8|j29[| , or by solar neutrino mixing at and beyond "maximal" values (uj > 7r/4) 



are currently being revisited (see [|30|-|3"2| and p3|-|35|, respectively) and demand an 
updated discussion. In particular, the case of maximal solar neutrino mixing (uj = 7r/4) 
appears now more "natural" and interesting in light of the Super-Kamiokande results on 
atmospheric neutrinos, which favor maximal <-> v T mixing (-0 = 7r/4).0 Therefore, it 
is important to go beyond the usual analyses in terms of the solar v mixing parameter 
sin 2 2uj (equivalent to take uj < 7r/4), which miss a potentially interesting region of the 
MSW solutions at uj > tt/4, as we shall see later (see also || for earlier discussions). 

Concerning the parameter 0, we do not include a priori the constraints coming from 
the CHOOZ reactor experiment f39j, which, together with the atmospheric neutrino data 



imply that sin 2 < few % [0)0 (^ ne exac ^ upper limit depending on the value of 



m 



on the confidence level chosen, and on the number of degrees of freedom in the oscillation 
model). In fact, we think that it is instructive to study the constraints on coming from 
solar neutrino data alone, as it has been done similarly for atmospheric neutrino data |f40| . 

In any case, observable effects on the MSW solutions can be generated even by small 
(few %) values of sin 2 0, which induce a fractional suppression in the solar neutrino survival 
probability P ee approximately equal to 2 x sin 2 0. In fact, if N e is the electron density profile, 
the 3v and 2v expressions for P ee are related by the expression 



Pl:{bm\uj, 



sin 



+ cos 4 0-P 2 ( T(5m 2 ,^ 



V e ^COS 2 <f)N e 



which, for small values of 0, gives roughly 



P; 



3v 



p. 



2v 



p2u 



-2 sin' 



(2) 



(3) 



Even within the stringent CHOOZ limits on sin 2 0, such variations of P ee can become as 
large as the uncertainties on the neutrino event rates, and thus cannot be neglected. 

The above motivations warrant the present work, which is organized as follows. In Sec. [II 
we present the theoretical and experimental ingredients of the analysis. In Sec. |ITT|we update 
the familiar two-family oscillation fit (0 = 0), and discuss in detail the features of the 2v 
MSW solutions. In Sec. we extend the analysis to three-flavor oscillations, and discuss 
how the 3z/ MSW solutions change for increasing values of 0. We summarize our results in 
Sec. |V[ The details of the statistical analysis are given in the Appendix. 



2 Solar v MSW oscillations at sin 2 2oj = 1 are nontrivial when Earth matter effects are included 
§j3f-|5|. This fact was not considered in Refs. |3q-p8|, where only matter effects in the Sun were 
studied. 
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II. STANDARD PREDICTIONS AND EXPERIMENTAL DATA 



In the present analysis, we use the so-called BP98 standard solar model [23] for the 
electron density in the Sun and for the input neutrino parameters (z/ e fluxes, spectra, and 
production regions), and compare the predictions to the experimental data for the following 
observables: total neutrino event rates, SK energy spectrum, and SK day-night asymmetry. 

Table | shows the latest results (and standard expectations) for the total neutrino event 
rates measured at Homestake ||19|| , Kamiokande p0| , SAGE |2_T|, GALLEX [22], and Super- 
Kamiokande (825 live days) 



12,13 



Since the SAGE and GALLEX detectors measure 
exactly the same quantity (up to a negligible difference in the detector latitude), their results 
are combined in a single (Ga) rate of 72.5±5.6 SNU. On the other hand, the Kamiokande and 
Super-Kamiokande data are treated separately (rather than combined in a single datum), 
since the two experiments, although based on the same v-e scattering detection technique, 
have rather different energy thresholds and resolution functions.^ 

The SK electron recoil energy spectrum and its statistical and systematic uncertainties 
(825 lifetime day, E e > 5.5 MeV) are graphically reduced from the 18-bin histograms shown 
by SK members in recent Summer '99 conferences P JT2]JT3[] . The corresponding numerical 
values have been already reported in detail in Ref. and are not repeated here. Our 
theoretical calculation of the binned spectrum properly takes into account energy threshold 
and resolution effects (see, e.g., the Appendix of Ref. |41[]). Standard 8 B |24| and hep f2~3 



neutrino spectra and fluxes are used, unless otherwise noted. Concerning the SK day-night 



asymmetry of the event rates, we use the latest measurement [13|: 



N-D 
N + D 



0.065 ±0.031 ±0.013 



(4) 



In the presence of 2u or 3z/ oscillations, the MSW effect in the Sun is computed as 
described in Ref. ||. The additional Earth matter effects are treated as in Ref. [27]. The 



X analysis of the theoretical and experimental uncertainties basically follows the approach 



developed in |42]] , with the necessary updates to take into account the BP98 SSM predictions 
and the energy spectrum information. Technical details about error estimates are given in 
the Appendix. 

We conclude this section by comparing the standard (no oscillation) predictions with 
the experimental data for the CI, Ga, and SK total rates. Figure [L] shows the 99% C.L. 
error ellipses for data and expectations in the planes charted by the (CI, Ga), (SK, Ga), and 
(SK, CI) total rates. The distance between observations and standard predictions makes 
the solar neutrino problem(s) evident. At present, such information is the main evidence for 
solar neutrino physics beyond the standard electroweak model; however, since the theoretical 



3 The inclusion of the Kamiokande rate is currently not decisive in shaping the MSW solutions, 
the SK rate being much more accurately measured. However, its addition to the CI, Ga, and SK 
rates avoids a situation of "zero degrees of freedom" in the x 2 fit to the total rates for the 3u MSW 
case [three data (CI, Ga, SK) minus three free parameters (5m 2 ,cu, cf))]. In such a situation, the 
value of Xmin would not have a well-defined likelihood. 
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errors are dominant — as far as total rates are concerned — no substantial improvements can 
be expected by a reduction of the experimental errors. Conversely, decisive information is 
expected from the SK spectrum and day-night asymmetry, but no convincing deviation has 
emerged from such data yet. Therefore, it is not surprising that, in oscillation fits, the total 
rates mainly determine allowed regions, while the SK spectrum and day-night asymmetry 
determine excluded regions. 



III. TWO-FLAVOR MSW OSCILLATIONS 

Figure |2| shows the results of our 2v MSW analysis of the data discussed in the previous 
section, shown as confidence level contours in the (6m 2 , sin 2 2cj/ cos2u;) plane. The choice 
of the variable sin 2 2ui / cos 2a;, rather than the usual sin 2 2uj, allows an expanded view of the 
large mixing region. In each of the six panels, we determine the absolute minimum of the x 2 
and then plot the iso-x 2 contours at x 2 ~ Xmin = 4.61, 5.99, and 9.21, corresponding to 90%, 
95%, and 99% C.L. for two degrees of freedom (the oscillation parameters). In fits including 
the total rates, there is a global x 2 minimum and two local mimima; such minima, and the 
surrounding favored regions, are usually indicated as MSW solutions at small mixing angle 
(SMA), large mixing angle (LMA), and low Sm 2 (LOW). Additional information on such 
solutions is reported in Tab. JTj]. 

Concerning the statistical interpretation of the x 2 values, a remark is in order. One 
can attach confidence levels to x 2 values in two different ways, depending on the choice 



between hypotheses tests and parameter estimation J43|. If one is interested in testing the 
goodness of the MSW hypothesis a priori, then one should compare the absolute Xmm with 
a number of degrees of freedom (A^f) calculated as number of data minus number of free 
parameters. The corresponding probability P is given in the last column of Tab. [TT] where, 
for completeness, P is reported also for the other two local minima and for the no oscillation 
case. If the MSW hypothesis is accepted, then the MSW parameter estimation involves only 
X 2 differences with respect to the global minimum, and the appropriate value of A DF to use 
is the number of free parameters in the model (Adf = 2), as anticipated for Fig. |2|. 

The first panel of Fig. |2| refers to the fit to the total rates only. The three x 2 minima are 
indicated by dots. The absolute minimum is reached within the SMA solution (Xmin = 1-08), 
which represents a very good fit to the data. The LMA solution is also acceptable, while 
the LOW solution gives a marginal fit (see also the upper three rows of Tab. [1]]). The SK 
data on the day-night asymmetry (second panel) and energy spectrum (third panel) exclude 
large regions in the mass-mixing parameter space; but are unable to (dis)prove any of the 
three solutions, which in fact are present also in the global fit to all data (fifth panel; see 
also the middle three rows of Tab. 0). 

The spectrum information is somewhat sensitive to the (uncertain) value of the hep 
neutrino flux; for instance, an enhancement by a factor 20 helps to fit the high-energy part 
of the SK spectrum ^0], and thus it produces a reduction of the excluded regions in the 
mass- mixing plane (fourth panel in Fig. 0), and a corresponding slight enlargement of the 
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globally allowed regions (sixth panel; see also the lower three rows of Tab. 

From the results of Fig. |2| and of Tab. |J it appears that the inclusion of the day-night and 
spectral information can significantly change the C.L.'s associated to each of three solutions 
SMA, LMA, and LOW. In order to understand better the role of difference pieces of data, 
we show in Fig. ^| the comparison between data and predictions for the Xmm point of the 
total rate fit (SMA solution at best fit, first row of Tab. H), in the same coordinates as in 
Fig. [TJ. Analogously, Figs. [| and [5] show the analogous comparison for the LMA and LOW 
solutions (second and third row of Tab. U). Figure |] shows the spectral information for the 
three local minima in the global fit to all data with standard hep flux (upper panel) and 
with enhanced (20 x) hep flux (lower panel). 

From Figs. it appears (as far as total rates are concerned) that the SMA solution 
represents a very good fit to the SK, Ga, and CI data, that the LMA solution underestimates 
slightly the Ga and SK data, and that the LOW solution underestimates the Ga rate and 
overestimates the CI rate. This explains the ordering in the likelihood of such solutions in 
Tab. |5 P(SMA) > P(LMA) > P(LOW). On the other hand, Figure § shows that the 
bulk of the SK observed spectrum is basically consistent with a flat shape (except for the 
two highest-energy, low-statistics bins), and therefore tends to favor the LMA and LOW 
solutions, rather than the steadily increasing spectrum predicted by the SMA solution. 
This tendency is slightly enhanced in the case of large hep neutrino flux (lower panel of 
Fig. H) since in this case the LMA and LOW solutions can fit somewhat better the spectrum 
endpoint. In other words, the inclusion of the spectrum in the global analysis tends to 
compensate the different likelihoods of the LMA and LOW solutions with respect to the 
SMA solution. Moreover, the slight excess of observed nighttime events [Eq. (|])] adds extra 
likelihood to the LMA solution, so that in the global fit one has P(SMA) ~ P(LMA) > 
P(LOW) (see Tab. [TT|), with the LMA solution even more likely than the SMA solution in 
the case of large hep flux. 

From the previous discussion, the following situation emerges for the three MSW solu- 
tions SMA, LMA, and LOW. None of them can be excluded at 99% C.L. by the present 
experimental data. Different pieces of data give indications that are not as consistent as it 
would be desirable: the total rate information favors the SMA solution, the spectral data 
favor the LMA and LOW solutions, and the day-night data favor the LMA solution. In a 
global fits, the three solutions have comparable likelihoods. Although such solutions are sub- 
ject to change shape and likelihood as more accurate experimental data become available, 
no dramatic improvement can be really expected in their selection, unless: 1) the theoretical 
uncertainties on the total rates are reduced to the size of the corresponding experimental 
uncertainties; 2) the total errors associated to the SK spectrum and day-night measurement 
are significantly reduced (by, say, a factor ~ 2); or 3) decisive results are found in new solar 
neutrino experiments such as the Sudbury Neutrino Observatory (SNO) the Gal- 



4 The SK Collaboration [jl2| derives a preliminary upper limit of ~ 15 to the hep/SSM flux ratio 
(in the absence of oscillations). This limit is weakened to ^ 30 in the presence of an oscillation 
effect of ~ 50%. Therefore, we can assume 20 x hep as a relevant, representative case of large hep 
flux. 
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lium Neutrino Observatory p6[p[, BOREXINO [|8],[|7|, an d KamLand (49|. Any of these 



conditions require a time scale of a few years at least; the same time scale should then be 
expected in order to (patiently) single out one of the three MSW solutions (SMA, LMA, or 
LOW). 

Another aspect of the LMA and LOW solutions emerging from Fig. ^| is their extension 
to large values of the mixing angle (sin 2 2uj — > 1), which are often assumed to be realized only 
through the vacuum oscillation solutions. Since the possibility of nearly maximal (^1,^2) 
mixing for solar neutrinos has gained momentum after the SK evidence for maximal (v^, v T ) 
mixing (sin 2 2ip ~ 1), it is interesting to study it in detail by dropping the usual "2a;" 
variable and by exploring the full range uj G [0, 7r/2], as it was done earlier in |J. The 
subcase uj = ir/A will receive special attention in the next section. 



IV. THREE-FLAVOR MSW OSCILLATIONS 

As stated in Sec. |, for large values of m 2 (^> 10~ 4 eV 2 ) the parameter space relevant for 
?>v solar neutrino oscillations is spanned by the variables (5m 2 , uj, 0). As far as uj is taken in 
its full range [0, 7r/2], one can assume 5m 2 > 0, since the MSW physics is invariant under 
the substitution (5m 2 , uj) — > (—5m 2 ,7i/2 — uj) at any 0.Q 

For graphical representations, we prefer to use the mixing variables (tan 2 uj, tan 2 0) intro- 
duced in [||], which properly chart both cases of small and large mixing. The case tan 2 = 
corresponds to the familiar two-family oscillation scenario, except that now we consider 
also the usually neglected case uj > 7r/4 (tan 2 u; > 1). For each set of observables (rates, 
spectrum, day-night difference, and combined data) we compute the corresponding MSW 
predictions and their uncertainties, identify the absolute minimum of the y 2 function, and 
determine the surfaces at x 2 — Xmin = 6.25, 7.82, and 11.36, which define the volumes con- 
straining the (5m 2 , tan 2 uj, tan 2 0) parameter space at 90%, 95%, and 99% C.L. Such volumes 
are graphically presented in (5m 2 , tan 2 a;) slices for representative values of tan 2 0. 

Figure ^ shows the results for the fit to the total rates only. The absolute minimum 
(Xmm = 0.7) is reached within the SMA solution at nonzero tan 2 0, as reported in the first 
row of Tab. The preference for 7^ is, however, not statistically significant, since 
x 2 min ((j) = 0) - Xmin = 1-08 - 0.70 = 0.38 only. The qualitative behavior of the 3z/ MSW 
solutions in Fig. [7] is unchanged with respect to the earlier analysis in 0: for increasing 
tan 2 0, the LOW solution moves to slightly larger values of tan 2 a; and then disappears, 
while the SMA and LMA regions tend to merge in a single, broad solution which slowly 
disappears at large values of 0. In general, the MSW allowed region(s) become less structured 
for increasing 0, due to the the decreasing energy dependence of the v e survival probability 
||. In the first panel of Fig. [7] it is interesting to note that both the LMA and the LOW 
solutions are compatible at 99% C.L. with maximal mixing (tan 2 u; = 1). This compatibility 
persists for the LOW solution at small values of tan 2 0, while it is lost for the LMA solution, 
which migrates towards smaller mixing angles. As for the two-family case discussed in the 



5 Only for solar neutrino oscillations in vacuum, the further symmetry uj — > tt/2 — to (at any 5m 2 
and E [0, n/2]) allows to take dm 2 > with oj in the restricted range [0, 7r/4]. 
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previous section, the addition of the spectrum information (and, to a lesser extent, of the 
day-night asymmetry) is expected to produce significant changes in the C.L.'s associated to 
the MSW solutions in Fig. 0. 

Figure || shows the region excluded by the SK day-night asymmetry, in the same repre- 
sentation as in Fig. [7]. Such region is generally far from the allowed regions in Fig. [7], and thus 
the day-night asymmetry contributes only marginally to shape the MSW solutions, except 
at low values of tan 2 0, where it cuts the lower part of the LMA solution and disfavors the 
rightmost part of the SMA solution. The region excluded by the SK day-night asymmetry 
can extend beyond maximal mixing (tan 2 a; > 1), as it was observed earlier in Ref. [Q for 
the corresponding Kamiokande variable. 

Figure || shows instead the regions excluded by the SK spectrum, which plays a relevant 
role in disfavoring the zone where the SMA and LMA solutions tend to merge (see the 
six middle panels of Fig. |^). In fact, for 5m 2 ~ 10~ 4 eV 2 the normalized SK spectrum 
is predicted to have a large negative slope (see, e.g., |P7||), contrary to the experimental 
results. Relatively large positive values for the spectrum slope are excluded in the left 
diagonal band shown in the upper panels of Fig. |9[. Both positive and negative values for 
the spectrum slope can instead occur in the excluded region where the Earth regeneration 



effects are relevant |27|| . For increasing 0, the spectral information becomes less effective 
in disfavoring zones of the parameter space, since the energy-dependence of the electron 
survival probability becomes weaker. Therefore, only the strongest distortion effect (at 
Sm 2 ~ 10 -4 eV 2 ) survives at large tan 2 0. Finally, as in the 2u case, the excluded regions 
are slightly shrinked if large hep flux values are assumed (not shown). 

Figure [10] shows the combined fit to all data, which is one one of the main results of this 
work. As for the total rates, the minimum \ 2 (second row of Tab. pl|) is reached within the 
SMA solution and shows a very weak preference for nonzero values of (tan 2 ~ 0.1). By 
comparing Fig. [IFJ with Fig. [?], it can be seen that the SK spectrum excludes a significant 
fraction of the solutions at 5m 2 ~ 10 -4 eV 2 , including the upper part of the LMA solution 
at small 0, and the merging with the SMA solution at large 0. In particular, at tan 2 = 0.1 
the 95% C.L. upper limit on 5m 2 drops from 2 x 10~ 4 eV 2 (rates only) to 8 x 10~ 5 eV 2 (all 
data). This indication tends to disfavor neutrino searches of CP violation effects, since such 
effects decrease with Sm 2 /m 2 at nonzero values of 0. 

Figure |1(] also shows that the inclusion of the present spectrum and day-night informa- 
tion, rather than helping in constraining tan 2 0, tends to slightly weaken its upper bound: 
for instance, the region allowed at 95% C.L. in the last panel of Fig. [TO is larger than the 
corresponding one in Fig. 0. In fact, the information coming from total rates prefers some 
energy dependence in the survival probability P ee , and thus disfavors large values of [see 
Eq. (0)]. However, since no clear indication for dP ee /dE v ^ emerges from the bulk of the 
SK spectrum, its addition in the fit makes large values of slightly more "acceptable." In 
addition, at tan 2 0^1, the lower part of the broad region allowed by total rates (last three 
panels of Fig. ^ becomes more consistent with the (weak) indication for a nonzero night-day 
difference [Eq. (0)], thus giving some extra allowance for large values of 0, as shown in the 



last three panels of Fig. [10 



Figure |Tj] is analogous to Fig. ID, the only difference being the flux of hep neutrinos 
(increased by a factor of 20 with respect to BP98). Since a large hep flux allows a better fit 
to the SK spectrum endpoint at any value of 0, the global fit is somewhat improved (third 
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row of Tab. |TTT[) , and the MSW solutions are slightly more extended than in Fig. [10]. As 



for the 2v case, the best-fit point migrates from the SMA to the LMA solution. In general, 
the comparison of Figs. [10] and [11] shows that the (uncertain) value of the hep neutrino flux 
is not crucial in shaping the current 3i> MSW solutions (although it might become decisive 
with more accurate spectral data). 

Figure [12] shows the values of \ 2 ~ Xmm as a function of tan 2 <fi for unconstrained Sm 2 
and uj, coming from the fit to total rates and to all data. The values of Xmm can be read 
from Tab. ITJ. This figure displays more clearly some anticipated features of the current Zv 
MSW solutions (for standard hep flux): 1) a slight preference for nonzero values of tan 2 0; 
2) the existence of an upper bound on 0, slightly stronger in the fit to total rates (0 < 55° 
at 95% C.L.) than in the global fit (0 ^ 58°). In the case of enhanced (20 x) hep flux we 
find that: 1) the absolute minimum of \ 2 is reached at tan 2 = within the LMA solution; 
2) for fixed 0^0, the local minimum migrates to the SMA solution when tan 2 > 0.05; 
and 3) the 95% C.L. upper limit on tan 2 is slightly weakened to ~ 59°. 

The 95% C.L. upper bound on coming from solar neutrino data alone (0 < 55°-59°) 
is consistent with the one coming from atmospheric neutrino data alone (0 ^ 45°) ||40|| , 
as well as with the upper limit coming from the combination of CHOOZ and atmospheric 
data (0 ^ 15°) |35|||(J. This indication supports the possibility that solar, atmospheric, and 
CHOOZ data can be interpreted in a single three-flavor oscillation framework [p^^Of. In this 
case, the CHOOZ constraints on exclude a large part of the 3z/ MSW parameter space 
(basically all but the first two panels in Figs. ]^-[T0|). 

However, even small values of can be interesting for solar v phenomenology, as an- 
ticipated in the comment to Eq. (|3]). Figure ^ shows, for instance, how the partial and 
global fits to neutrino event rates (CI, Ga, and SK experiments) are modified when passing 
from tan 2 = to tan 2 = 0.06. Each of the three experiments, taken separately, allows 
uj < 7r/4. In the combination, the LOW and LMA solutions "touch" the u = tt/4 line at 
tan0 = 0, while for nonzero only the LOW solution remains consistent with u = tt/A 
(and, to some extent, with u > 7r/4). This fraction of the LOW solution, usually ignored in 
analyses using the sin 2 2uj variable, can be relevant for model building. In particular, there 
is currently great interest in models predicting exact or nearly bimaximal mixing, namely, 
(uj, ip, 0) ~ (tt/4, 7r/4, 0) (see, e.g., |50| , |5T| and refs. therein). Such models imply, for the neu- 
trino mixing matrix U a i, that U 2 X ~ U 2 2 and that U 2 3 ~ 0. Recent theoretical predictions 
for (small) nonzero values of U 2 3 are discussed, e.g., in |52 . 

Figure [T3] shows the section of the volume allowed in the 3u MSW parameter space, for 
uj = 7r/4 (maximal mixing), in the mass-mixing plane (Sm 2 , sin 2 0). All data are included 
in the fit (with standard hep flux. It can be seen that both the LMA and LOW solutions 
are consistent with maximal mixing (at 99% C.L.) for sin 2 = U 2 3 = 0. Moreever, the 
consistency of the LOW solution with maximal mixing improves significantly for U 2 3 ~ 0.1, 
while the opposite happens for the LMA solution. This gives the possibility to obtain nearly 
bimaximal mixing (uj = ip = it/A with small) within the LOW solution to the solar 
neutrino problem — an interesting possibility for models predicting large mixing angles. 
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V. SUMMARY 



We have performed a thorough analysis of the 2u and 3z/ MSW solutions of the solar 
neutrino problem, discussing the information coming from the total rates, as well as from 
the SK energy spectrum and day-night asymmetry. The global fit to the data puts an 
upper bound to the #13 mixing angle, consistent with atmospheric and reactor oscillation 
searches. A fraction of the MSW solutions extends at and beyond maximal (ui, U2) mixing 
(612 > 7r/4), especially for the so-called LOW solution. Therefore, a novel realization of 
(nearly) bimaximal mixing models appears possible. 
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APPENDIX: STATISTICAL ANALYSIS 



In this appendix we present an updated version of the x 2 statistical analysis of neutrino 



event rates discussed in Ref . |42| , which has been widely used in many other analyses of the 
solar neutrino problem, including the one performed by the Super-Kamiokande Collabora- 
tion in In particular, we take into account the most recent estimates for the relevant 
BP98 standard solar model ingredients and for their uncertainties 



23jj. Unless otherwise 



noted, the notation is the same as in Ref. 
details. 



to which the reader is referred for further 



as calculated in the BP98 Standard Solar 



Table [TV] reports the solar neutrino fluxes 
Model The corresponding energy-averaged interaction cross sections Cy for the j- 

th detector are given in Tab. |V|, taking into account updated energy spectra ||23|.|24|] and 
cross sections p^j26|,|25[] . The associated la relative uncertainties AlnCy are reported in 



Tab. [Vl|. Such tables allow the calculation of the theoretical uncertainties related to the 
detection process. 

The uncertainties of the SSM v fluxes are mainly related on eleven basic ingredients 
Xfc: five astrophysical S'-factors, the Sun luminosity, the metallicity Z/X, the Sun age, the 
opacity, the element diffusion, and the 7 Be electron capture rate Cb c - The last two sources 
of uncertainties represent new entries with respect to the list given in ||42|| . The "diffusion 
error" is estimated by taking differences between quantities calculated with and without 
diffusion f23[. The uncertainty in the 7 Be(e~, z/ e ) 7 Li capture rate [^3|] affects (inversely) the 



8 B flux by changing the 7 Be density for the (much slower) competing reaction 7 Be(p, 7) 8 B. 
The fractional la uncertainties of the eleven input ingredients (AlnXfc) are reported in 
Tab. |VTJ 

Table |VIII| reports the matrix of logarithmic derivatives <91n0j/<91n Xk, which pa- 
rameterize the SSM "response" to small variations in the input ingredients, and which are 
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crucial to estimate correctly the correlations among the neutrino flux or event rate uncer- 



tainties p2| . As remarked in ||2| for the opacity uncertainty, we have split the known flux 
errors Aln0j = a^AlnX^ due to opacity, diffusion, and 7 Be electron capture, into two 
arbitrary factors, a ik and AlnX fc (k = 9, 10, 11), so as to treat homogeneously all sources 
of uncertainties in the \ 2 statistics. 

Given the previous ingredients, the standard event rates Rj are calculated as: 



and the associated theoretical error matrix crj^ reads: 



2 



(Al) 



(A2) 



The off-diagonal entries of erj • induce strong correlations among the theoretical uncertain- 
ties for different experiments, as evident in Fig. [l|. In the presence of oscillations (Figs. |3|-[5|), 
the standard values of the partial rates are simply replaced by the corresponding oscil- 
lated values. 

We remind that the entries of the matrix are not independent, being constrained by 



the luminosity sum rule [42|. At equilibrium, the standard neutrino fluxes must obey the 
energy conservation relation ffl: 



£(Q/2 - Et) <p t = K 



(A3) 



where is the average neutrino energy for the z-th neutrino flux, K = 0.853 x 10 12 
MeV cm -2 s _1 is the solar luminosity, and Q = 26.73 MeV is the overall Q-value for the 
fusion reaction 4p — > a + 2u e + 2e + (independent from the intermediate reaction chain). 

Non-equilibrium corrections, which are especially relevant for CNO fluxes, modify this 
relation into [l54l|55l 



K 



(A4) 



where the £j values are given in Tab. |TX|. By taking partial derivatives of the above relation 

(A5) 



with respect to the input ingredients X^, one gets the sum rule | 42| : 

J2i & <t>i Oiik 



6, 



A-6 



which is satisfied with good approximation (~ 2% level). 

Similarly, if the flux variations are approximately parametrized as power laws of the 
central central temperature T of the Sun, 



rp \f3i 

AH = 0,(1,) ( -J 



(A6) 



the temperature exponents (3i become constrained by the sum rule fl4*2",5"l| 
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^iMl = (A7) 



which, using the exponents in Ref. ||54|| , is satisfied only at the ~ 10% level. However, such 
exponents can be easily adjusted within their uncertainties — > /%) so as to satisfy exactly 
the above sum rule. Sets of temperature exponents are given in Tab. |X|. The set of /3^s is 
particularly relevant for analyses of the solar neutrino problem with "unconstrained" central 
temperature (not performed in this work). 

Since the above tables and equations refer only to the analysis of the total neutrino event 
rates, a final remark is in order about the fit to the SK energy-angle spectral information. 
For MSW analyses, the maximum amount of information about the SK solar neutrino signal 
is contained in the double-differential distribution of events in nadir angle (77) and electron 
recoil energy (E e ). However, the (77, E e ) distribution is not yet available outside the SK 
collaboration, which has released only the two projected distributions, namely, the (angle- 
averaged) energy spectrum and the (energy-averaged) nadir distribution. The experimental 
errors of such projections (including statistical fluctuations) are a priori strongly correlated, 
but the published information is not sufficient to recover such correlations. For such reasons, 
while we include the complete information in E e above 5.5 MeV (given as an 18-bin spectrum 
in 0,[E||T3|]), we prefer to use only a minimal (but relevant) information in 77, namely, the 
night-day asymmetry (N — D)/ (N + D). 

The x 2 function for the energy spectrum is then constructed in the same way as in 
Ref. |17]], namely, by separating uncorrelated and correlated error components in the 18 x 18 



error matrix, with off-diagonal correlations equal to +1. A numerical table of errors is 



explicitly reported in [I7| for the 825 day SK spectrum. Double counting of the total rate 
information in SK is avoided by minimizing the x 2 with respect to the overall spectrum 



normalization, taken as a free parameter ||17|| . In the global fit, we then sum up the x °f 



the spectrum, the \ 2 of the total rates, and the (trivial) x 2 °f the day-night asymmetry 
measurement. A similar construction of the global x 2 function is also adopted in WB\. 
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TABLES 



TABLE I. Solar neutrino event rates observed in the five experiments Homestake, Kamiokande, 
SAGE, GALLEX, and Super-Kamiokande (825 days, E e > 6.5 MeV), as compared with the the- 
oretical expectations from the BP'98 standard solar model |23|]. The quoted uncertainties are at 
lcr. In our analysis, the GALLEX and SAGE experimental rates are combined in quadrature in a 
single (Gallium) rate of 72.5 ± 5.6 SNU. 



Experiment 



Experimental rate 



Theoretical rate 



fir 

-1.0 
+1.0 
-0.7 

+8 



Units 



Ref. 



Homestake 

Kamiokande 

SAGE 

GALLEX 

Super-Kamiokande 



2.56 ±0.16 ±0.16 
2.80 ±0.19 ±0.33 



nn 9 +7.2 +3.5 
°' - z -7.0 -3.0 

77.5 ± 6.2 It? 



2.45 ±0.04 ±0.07 



7.7 
5.15 
129 



129 jj 
°- i0 -0.7 



SNU 
10 6 cm" 2 s" 1 

SNU 

SNU 
10 6 cm" 2 s- 1 



IS] 

m 
m 

Tl 



TABLE II. 2v MSW analysis: Absolute x 2 values in selected points of the parameter space 
(5m 2 , sin 2 2cj/ cos 2ui), corresponding to local x 2 minima (SMA, LMA, LOW solutions) and to the 
no-oscillation case. Upper three rows: fit to total rates only (Cl±Ga±K±SK data). Middle three 
rows: fit to total rates (4 data), day-night asymmetry (1 datum), and SK energy spectrum (18 
data with one adjustable normalization parameter). Lower three rows: as for the middle rows, but 
with hep flux increased by a factor 20. The number of degrees of freedom A^df is equal to the 
number of data (either 4 or 4 + 1 + 18 — 1) minus 2 (the oscillation parameters). The last column 
corresponds to the probability of having a worse x 2 fit a priori. 



Observables 


Solution 


5m 2 (eV 2 


) 


sin 2 2lo/ cos 2uj 


x 2 


^(x'^df) 




SMA 


5.89 x 1(T 


-6 


5.62 x 10~ 3 


1.08 


58.3 % 


Rates only 


LMA 


1.91 x 10" 


-5 


1.62 


4.59 


10.1 % 


(AT DF = 2) 


LOW 


1.07 x 10" 


-7 


3.48 


8.24 


1.6 % 




no osc. 










60.1 


9 x 10" 14 


Rates, 


SMA 


5.37 x 10" 


-6 


6.50 x 10" 3 


27.3 


12.7 % 


N- D/N + D, 


LMA 


2.55 x 10- 


-5 


1.62 


27.8 


11.4 % 


spectrum 


LOW 


1.07 x 10- 


-7 


2.87 


30.8 


5.8 % 


(N DF = 20) 


no osc. 










85.8 


4 x lO" 10 


Rates, 


SMA 


5.20 x 10- 


-6 


6.20 x 10" 3 


25.6 


17.9 % 


N- D/N + D, 


LMA 


2.94 x 10- 


-5 


1.96 


24.6 


21.7 % 


spectrum [hepx 20] 


LOW 


1.24 x 10- 


-7 


2.87 


27.8 


11.4 % 


{N DF = 20) 


no osc. 










71.9 


9 x 10~ 8 
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TABLE III. 3v MSW analysis: Absolute minima of the function ^(Sm 2 , tan 2 u>, tan 2 <j>) and 
their associated probability P. The number of degrees of freedom is iVrjF = 1 for the fit to total 
rates only and JVdf = 19 for the global fit. 



Observables 


5m 


2 (eV 2 ) 


tan 2 u> 


tan 2 cf) 


x 2 - 

A-min 


P(X 2 ,N DF ) 


Rates only 


9.8 


x 10' 6 


8.1 x 10" 4 


0.1 


0.70 


40.3 % 


All data 


9.8 


x 10~ 6 


7.1 x 10~ 4 


0.1 


27.0 


10.5 % 


All data [hep x 20] 


2.8 


x 10~ 5 


0.37 





24.5 


17.8 % 



TABLE IV. Neutrino fluxes <fc (cm" 2 s" 1 ) from the BP98 SSM ||. 



pp 


pep 


hep 


Be 


B 


N 


O 


F 


xlO 10 


xlO 8 


xlO 3 


xlO 9 


xlO 6 


xlO 8 


xlO 8 


xlO 6 



5.94 1.39 2.10 4.80 5.15 6.05 5.32 6.33 



TABLE V. Energy-averaged cross sections Cij. Units are such that the expected partial rates 
Rij = Cij4>i are directly given in SNU for the Ga and CI experiments, and as a ratio to the total 
BP98 predicted rate for Super-Kamiokande {E e > 6.5 MeV) and Kamiokande. 





PP 


pep 


hep 


Be 


B 


N 


O 


F 




xlO" 9 


xlO" 8 


xlO" 6 


xlO" 9 


xlO" 6 


xl0~ 9 


xl0~ 8 


xl0~ 8 


Ga 


1.17 


2.04 


7.14 


7.17 


2.40 


6.04 


1.14 


1.14 


CI 





0.16 


4.26 


0.24 


1.14 


0.17 


0.07 


0.07 


K 








0.595 





0.194 











SK 








0.514 





0.194 












TABLE VI. la relative errors A In Cij of the energy averaged cross sections. The uncertainties 
of the (radiatively corrected) v-e scattering cross section in K and SK are negligible. 





PP 


pep 


hep 


Be 


B 


N 





F 


Ga 


0.023 


0.170 


0.320 


0.070 


0.320 


0.060 


0.120 


0.120 


CI 





0.020 


0.037 


0.020 


0.032 


0.020 


0.020 


0.020 


K, SK 



























TABLE VII. la relative errors A In of the relevant SSM input parameters. 



5*11 


S33 


5*34 


5*1,14 


Sn 


Lum 


Z/X 


Age 


Opa 


Diff 


Cbc 


0.017 


0.060 


0.094 


0.143 


0.106 


0.004 


0.033 


0.004 


0.02 


0.02 


0.02 



14 



TABLE VIII. Logarithmic derivatives a,ik = dlnfa/dlnX^, parameterizing the response of 
neutrino fluxes fa to variations in the SSM input parameters Xf;. 





pp 


pep 


hep 


Be 


B 


N 


O 


F 


5*11 


+0.14 


-0.17 


-0.08 


-0.97 


-2.59 


-2.53 


-2.93 


-2.94 


-533 


+0.03 


+0.05 


-0.45 


-0.43 


-0.40 


+0.02 


+0.02 


+0.02 




-0.06 


-0.09 


-0.08 


+0.86 


+0.81 


-0.05 


-0.05 


-0.05 


Si, 14 


-0.02 


-0.02 


-0.01 


+0.00 


+0.01 


+0.85 


+1.00 


+0.01 


Sn 


+0.00 


+0.00 


+0.00 


+0.00 


+1.00 


+0.00 


+0.00 


+0.00 


Lum 


+0.73 


+0.87 


+0.12 


+3.40 


+6.76 


+5.16 


+5.94 


+6.25 


Z/X 


-0.08 


-0.17 


-0.22 


+0.58 


+1.27 


+1.86 


+2.03 


+2.09 


Age 


-0.07 


+0.00 


-0.11 


+0.69 


+1.28 


+1.01 


+1.27 


+1.29 


Opa 


+0.14 


+0.24 


+0.54 


-1.38 


-2.62 


-1.67 


-2.05 


O 1 Q 

— z.lo 


Diff 


+0.13 


+0.22 


+0.13 


-0.90 


-2.00 


-2.56 


-2.75 


-2.75 


CBe 


+0.00 


+0.00 


+0.00 


+0.00 


-1.00 


+0.00 


+0.00 


+0.00 


TABLE IX. Neutrino enerj 


ry parameters (in MeV) relevant for the luminosity 


sum rules 


[Eqs. (Xjj-55)]. 


















PP 


pep 


hep 


Be 


B 


N 


O 


F 




0.261 


1.442 


9.625 


0.813 


6.735 


0.706 


0.994 


1.000 


Q/2- 


Ei 13.10 


11.92 


3.74 


12.55 


6.63 


12.66 


12.37 


12.37 




13.10 


11.92 


10.17 


12.55 


6.63 


3.457 


21.57 


2.36 



TABLE X. Temperature exponents Pi and their ±lcr errors (from p4],p|). The exponents Pi 
represent our best-fit adjustments, which satisfy exactly the luminosity sum rule [Eq. ( |A7[ )]. 



PP 


pep 


hep 


Be 


B 


N 


O 


F 


Pi -1.1 ±0.1 


-2.4 ±0.9 


4.5 ± 1.5 


10 ±2 


24 ±5 


24.4 ±0.2 


27.1 ±0.1 


27.8 ±0.1 


Pi -1.1039 


-2.407 


4.5 


8.79 


23.996 


24.4 


27.099 


27.8 
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FIGURES 

FIG. 1. The solar neutrino deficit, shown as a discrepancy between data and expectations in 
the gallium (Ga), chlorine (CI), and Super-Kamiokande total event rates. In each plane, the error 
ellipses represent 99% C.L. contours for two degrees of freedom (i.e., A% 2 = 9.21). The projection 
of an ellipse onto one of the axis gives approximately the ±3<7 range for the corresponding rate. 
Data and expectations refer to Table |. The correlation of SSM errors is calculated as in the 
Appendix. 

FIG. 2. Two-generation MSW solutions to the solar neutrino problem. The upper four 
panels correspond to the following separate fits to data subsets: total rates (Cl+Ga+K+SK); Su- 
per-Kamiokande night-day asymmetry N — D/N + D; Super-Kamiokande electron energy spectrum 
with standard hep neutrino flux; Super-Kamiokande spectrum with enhanced (20 x) hep neutrino 
flux. The two lower panels show the results of global fits to all data. The thin solid, thick solid, 
and dashed curves correspond to x 2 ~ Xmm = 4.61, 5.99, and 9.21. The positions of the local x 2 
minima in fits including the total rates are indicated by dots. See also Tab. [ij. 

FIG. 3. The SMA solution at best fit (total rates only, first row of Tab. |0), compared with the 
experimental data, in the same planes as in Fig. |lj. Note the excellent agreement between theory 
and observations. 

FIG. 4. As in Fig. || but for the LMA solution at best fit (total rates only, second row of 
Tab. 0). 

FIG. 5. As in Fig. ||, but for the LOW solution at best fit (total rates only, third row of 
Tab. 0). 



FIG. 6. Comparison of data and predictions for the Super-Kamiokande spectrum shape. The 
theoretical spectra normalization is taken free in the fit. The upper (lower) panel corresponds to 
the case of standard (20x) hep flux. The SMA, LMA, and LOW spectra are calculated in the 
global best fit points reported in Tab. || (middle rows for the upper panel and lower rows for the 



lower panel). The SK data are reported from PJl2|,13|, and the error bars refer to statistical and 
total experimental errors. 



FIG. 7. Three-flavor MSW oscillations: global fit to Cl+Ga+K+SK rates in the 
(5m 2 , tan 2 u, tan 2 <p) parameter space. The favored regions in each panel correspond to sections of 
the volume allowed at 90%, 95%, and 99% C.L. ( X 2 - xLin = 6 - 25 > 7 - 82 > and 1L36 ) for represen- 
tative values of tan 2 cf>. 

FIG. 8. As in Fig. but for the fit to the Super-Kamiokande night-day asymmetry. The 
region inside the curves is excluded. 
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FIG. 9. As in Fig. [j], but for the fit to the Super-Kamiokande energy spectrum. The regions 
inside the curves are excluded. 



FIG. 10. Results of the global three-flavor MSW fit to all data. Notice that, in the first two 
panels, the 99% C.L. contours are compatible with maximal mixing (tan 2 u = 1) for both the LOW 
and the LMA solutions. 

FIG. 11. As in Fig. [lCl, but for the case of enhanced (20 x) hep flux. The allowed regions are 



slightly enlarged with respect to Fig. 10 



FIG. 12. Values of A% 2 as a function of 0, for unconstrained Sm 2 and tan 2 uj. At 95% C.L., 
the upper limit on (p is in the range 55°-59°, depending on the data used in the fit and on the 
value of the hep flux. 

FIG. 13. Regions allowed at 99% C.L. by the total rates only, for tan 2 = (solid curves) 
and tan 2 = 0.06 (dotted curves). For tan 2 = 0.06, the SMA and LMA solutions are slightly 
shifted to lower values of tan 2 to, while the LOW solution is shifted to higher values (including the 
value u> = tt /4). 

FIG. 14. Section of the allowed volume in the 3u parameter space in the plane (dm 2 , sin 2 0) , 
for the case of maximal (^1,^2) mixing (to = vr/4). For sin 2 = 0, both the LMA and LOW 
solutions are compatible with maximal mixing at 99% C.L. For small values of sin 2 0, the maximal 
mixing case favors the LOW solution. 
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Solar neutrino deficit, 1999 




FIG. |j. The solar neutrino deficit, shown as a discrepancy between data and expectations in 
the gallium (Ga), chlorine (CI), and Super-Kamiokande total event rates. In each plane, the error 
ellipses represent 99% CL. contours for two degrees of freedom (i.e., A% 2 = 9.21). The projection 
of an ellipse onto one of the axis gives approximately the =L3<7 range for the corresponding rate. 
Data and expectations refer to Table |. The correlation of SSM errors is calculated as in the 
Appendix. 
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2v MSW oscillations 
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FIG. y. Two-generation MSW solutions to the solar neutrino problem. The upper four panels 
correspond to the following separate fits to data subsets: total rates (Cl+Ga+K+SK); Super- 
Kamiokande night-day asymmetry N — D/N + D; Super-Kamiokande electron energy spectrum 
with standard hep neutrino flux; Super-Kamiokande spectrum with enhanced (20 x) hep neutrino 
flux. The two lower panels show the results of global fits to all data. The thin solid, thick solid, 
and dashed curves correspond to x 2 ~ Xmm = 4.61, 5.99, and 9.21. The positions of the local x 2 
minima in fits including the total rates are indicated by dots. See also Tab. ||. 



21 



MSW oscillations, SMA solution 
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FIG. |3[ The SMA solution at best fit (total rates only, first row of Tab. ||), compared with the 
experimental data, in the same planes as in Fig. |]. Note the excellent agreement between theory 
and observations. 
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MSW oscillations, LMA solution 
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FIG. |j. As in Fig. ||, but for the LMA solution at best fit (total rates only, second row of 
Tab. 0). 



23 



MSW oscillations, LOW solution 
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FIG. ^. As in Fig. [|], but for the LOW solution at best fit (total rates only, third row of Tab. U). 
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Super-Kamiokande electron energy spectrum 
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FIG. ^. Comparison of data and predictions for the Super-Kamiokande spectrum shape. The 
theoretical spectra normalization is taken free in the fit. The upper (lower) panel corresponds to 
the case of standard (20x) hep flux. The SMA, LMA, and LOW spectra are calculated in the 
global best fit points reported in Tab. || (middle rows for the upper panel and lower rows for the 
lower panel). The SK data are reported from PJl2],13|, and the error bars refer to statistical and 
total experimental errors. 
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3i/ MSW: Rates (CI + Ga + K + SK) 
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FIG. 0. Three-flavor MSW oscillations: global fit to Cl+Ga+K+SK rates in the 
(5m 2 , tan 2 u, tan 2 <fi) parameter space. The favored regions in each panel correspond to sections of 
the volume allowed at 90%, 95%, and 99% C.L. ( X 2 - Xmin = 6 - 25 > 7 - 82 > and 1L36 ) for represen- 
tative values of tan 2 <j>. 
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3i/MSW: N-D/N + D 
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FIG. ||. As in Fig. 0, but for the fit to the Super-Kamiokande night-day asymmetry. The region 
inside the curves is excluded. 
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3i/MSW: SK spectrum (hep/SSM=1) 
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FIG. ||. As in Fig. [?], but for the fit to the Super-Kamiokande energy spectrum. The regions 
inside the curves are excluded. 
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3i/MSW: Rates + N-D/N + D + Spectrum (hep/SSM = 1 ) 
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FIG. 10. Results of the global three-flavor MSW fit to all data. Notice that, in the first two 
panels, the 99% C.L. contours are compatible with maximal mixing (tan 2 u = 1) for both the LOW 
and the LMA solutions. 
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3i/ MSW: Rates + N-D/N + D + Spectrum (hep/SSM=20) 
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FIG. 11, As in Fig. 10, but for the case of enhanced (20 x) hep flux. The allowed regions are 



slightly enlarged with respect to Fig. 10 
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<P (deg) 




FIG. 12. Values of Ax 2 as a function of </>, for unconstrained 5m 2 and tan 2 to. At 95% C.L., the 
upper limit on (j) is in the range 55°-59°, depending on the data used in the fit and on the value of 
the hep flux. 
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FIG. |T^. Regions allowed at 99% C.L. by the total rates only, for tan 2 <ft = (solid curves) and 
tan 2 (j) = 0.06 (dotted curves). For tan 2 cj) = 0.06, the SMA and LMA solutions are slightly shifted 
to lower values of tan 2 u, while the LOW solution is shifted to higher values (including the value 

to = 7T/4). 
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FIG. 14. Section of the allowed volume in the 3v parameter space in the plane (8m 2 , sin 0), for 
the case of maximal (v\, v-i) mixing (cu = 7r/4). For sin 2 4> = 0, both the LMA and LOW solutions 
are compatible with maximal mixing at 99% C.L. For small values of sin 2 <ft, the maximal mixing 
case favors the LOW solution. 
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